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Abstract 

A new form to construct complex superpotentials that produce real energy spectra in 
supersymmetric quantum mechanics is presented. This is based on the relation between 
the nonlinear Ermakov equation and a second order differential equation of the Schrodinger 
type. The superpotentials so constructed are characterized by the Ermakov parameters in 
such a way that they are always complex-valued and lead to non-Hermitian Hamiltonians 
with real spectra, whose eigenfunctions form a bi-orthogonal system. As applications we 
present new complex supersymmetric partners of the free particle that are T^T-symmetric 
and can be either periodic or regular (of the Pbschl-Teller form). A new family of complex 
oscillators with real frequencies that have the energies of the harmonic oscillator plus an 
additional real eigenvalue is introduced. 


1 Introduction 

The factorization method is an ontstanding algebraic techniqne used to analyze the spec¬ 
tral properties of exactly solvable potentials in quantum mechanics HHa. Although the 
hrst record of the method can be traced back to the Fock article on configuration space 
and second quantization published in 1932 |1], and this was explicitly used by Dirac in the 
1935 edition of his book [5] “as a little stratagem to solve the spectral problem for the one¬ 
dimensional quantum oscillator” [3], the factorization is usually attributed to Schrodinger 
who wrote on the matter eight years after Fock [Hj. The main idea is to reduce the sec¬ 
ond order differential form of the one-dimensional Hamiltonian H = + V{x) to 

the product of two hrst order differential operators A, H, up to an additive constant e 
(proper units are assumed); that is, H = AB -\- e. In the harmonic oscillator case, the 
factors A and B are respectively the creation and annihilation a operators of the boson 
algebra [a, a^] = aa^ — a^a = 1, with e = ^ = i^o the ground state energy. Following 
Fock, H = N + ^ corresponds to the energy of a single mode of the electromagnetic 
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field and N = a^a accounts for the number of photons in that mode. The method is 
also in the kernel of supersymmetric quantum mechanics, introduced by Witten as a ‘toy 
model’ of what happens in quantum held theories [7] . This last connection arose renewed 
interest in the issues of factorization and facilitated the strengthening of supersymmetric 
quantum mechanics as an algebraic discipline that is interesting in itself [HHH] (for recent 
reviews see Nevertheless, the entire development of the factorization method 

“shows chronological gaps and inconsistencies; the ideas emerge, disappear and re-emerge 
again” [3]. 

The present work addresses the problem of exploring new prospects to the factorization 
method. In addition to the supersymmetric treatments already established to construct 
non-Hermitian Hamiltonians with either real eigenvalues or a combination of real and 
complex eigenvalues (see e.g. [T81 - I26] and references quoted therein), we use here the 
solutions of the Ermakov equation to construct complex-valued potentials that share their 
spectrum with a given solvable real potential. Common in references [IBH2S] is the idea 
of considering the eigenvalue equation Hu = eu, with u{x) not necessarily normalizable 
in Dom(l/) C R, as a seed repository of the Darboux transformations V = V — 2^(lnM). 
The latter allow the construction of new solvable potentials V in terms of the solutions 
of a given exactly solvable potential V and the ‘superpotentiah /3 = —^(Inw). Thus, 
V = V + 2^ can be constructed such that Dom(l/) = Dom(l/) and either cr(V) = cr(V) 
or (7(V) = o'(V) U {e}, with cr{A) the spectrum of A. The former case is referred to 
as ‘strictly isospectral’ and the second one as ‘almost isospectral’. In general, if both V 
and (3 are real-valued then the new potential V is real-valued. However, a large variety 
of Darboux deformations V can be achieved if the superpotentials f3 are complex-valued 
functions [I81 - I26] . Indeed, if V and {3 are respectively real and complex-valued functions 
then necessarily V is complex-valued. If both V and f3 are complex then V will be complex 
in general. 

To our knowledge, the current methods for constructing a complex-valued superpo¬ 
tential include at least one of the following ingredients: 

i) A normalizable complex function u belonging to the eigenvalue e G C of the non- 
Hermitian potential V E C. Usually Dom(U) = R"*" and m( 0) = 0 [H], but it can 
be extended to Dom(U) = R and the proper boundary conditions on u. 

ii) A complex linear combination f3 = f + ig oi real functions / and g that satisfy 
a set of four Riccati-like equations connecting the real and imaginary parts of the 
non-Hermitian potentials V and V [T9],|20]. 

hi) A complex linear combination u = ui + CU 2 , c G C, of the two fundamental eigen- 
solutions Ui,U 2 , associated to the eigenvalue e G R of the Hermitian potential U G R. 
Clearly Wim,U 2 ) = - ^U 2 ^ 0 [21]. 

iv) A complex analytic function u belonging to the eigenvalue e G C of the Hermitian 
potential U G R [221123] . 
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v) A Gamow-Siegert function u associated to the resonance e = E—^T of the Hermitian 
potential G G M [21H26]. 


However, the identihcation of exactly solvable complex potentials that can be generated 
by factorization (supersymmetric, Darboux or intertwining approaches) is, by far, not yet 
over. The purpose of the present work is to show that new classes of complex potentials 
with real spectrum are still achievable from the factorization method. With this aim we 
hrst review the connection between the solutions of the eigenvalue equatioiJi] Hu = eu 
and the solutions of the Ermakov equation 


(Pa 

dx^ 


Ma 


a-^ 


( 1 ) 


where Aq is a constant and M may be a function of x. This last nonlinear equation dates 
from a theorem proved by Ermakov in 1880 which states that if an integral of the equation 


fu 

dx"^ 


Mu 


( 2 ) 


is known, one can find an integral of the equation and vice versa, if a particular so¬ 
lution o/ ([I]) is known, we can find the complete solution of (^j [27], §20. Independent 
works by Pinney [2E], and Lewis and Riesenfeld [22], motivated the application of the 
Ermakov equation in physics [30H39] (for recent reviews see [101 ITT] ). Remarkably, the 
above mentioned relationship has been used to solve the quantum spectral problem of 
some potentials in terms of the so called Milne quantization condition [12TI11] . On the 
other hand, the relationship between ([1]) and the Riccati equation fulhlled by the real¬ 
valued superpotential jd = —^(ln(po), with (po the ground-state wavefunction of a given 
one-dimensional potential, has been indicated in [10] to connect the conventional super- 
symmetric quantum mechanics with the Ermakov approach (as discussed above, in this 
case all the Darboux deformations V = V — 2^(ln(po) are real-valued functions). In 
such a context, it has been shown that the Ermakov equation ([1]) obeys Darboux (su¬ 
persymmetric) covariance [16]. Here we consider stationary systems and show that, even 
for real potentials V(x) and real eigenvalues e, the solutions of the Ermakov equation 
m are useful to construct complex-valued solutions of Hu = eu. Using these last as 
‘transformation functions’ the superpotential (3 = — ^(Inn) can always be chosen to be 
complex-valued, so that all the new potentials V = V 2^ are complex. 

Our interest in complex-valued potentials includes the analogy between the spatially 
varying refractive index m{r) and the potential V (r) in the scattering of light and quantum 
particles respectively (in general, r is a vector in M^). Such analogy is usually drawn in 
terms of the identification np{f) = 1 — V{r)/[E], with [E] representing energy units [17] . 

^From now on we omit a global factor 1/2 in the Hamiltonian to adopt the simpler mathematical 
notation H = + V{x). The expressions for H and its eigenvalues as they usually appear in physics 

are easily recovered from our results by retrieving the 1/2 factor. For instance, in mathematical notation 
the potential and energy eigenvalues of the harmonic oscillator can be read as E (a;) = and En = 2n -|-1. 
The ‘physical expressions’ are respectively given by hM = and =n-\-\. 
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If the scattering medium is absorbing, the refractive index is complex m = n — ik, with 
nonzero imaginary part k > 0. Here n stands for the refractive index relative to the 
medium and k the absorption coefficient of the medium. We stress that k could be as small 
as to be omitted from calculations in a first approach, but this is never exactly zero for 
actual scattering media (see e.g. the experimental data reported in m for the values oik). 
Within this analogy, the non-Hermitian potentials V{f) = Vo(r() +iVi{f) appearing in the 
quantum “description of the scattering and the compound nucleus formation by nucleons 
impinging upon complex nuclei” [H] are referred to as ‘optical potentials’ [19] (see also 
‘atypical models’ in [3]). The simplest model considers Vo(r) as a radial rectangular well 
of depth U and width i?, and Vi(r) as a constant fraction of Vq, namely Vi = C,Vo with 
C > 0. Then Vo can be interpreted as the average potential in the nucleus and Vi as 
the parameter of absorption associated with the formation of the compound nucleus [H] 
(notice that Vi is the negative constant —^U in r < R). Using the factorization method 
with e = U —|r, T > 0, it has been shown that the well potential Vo(r) of [H] is indeed the 
seed of a family of optical potentials V{r) = Vo(r) + iVi(r), the imaginary part of which 
is not a constant but a function of the position that changes sign across different regions 
of [0,+oo) as r increases [2S]. Thus, Vi describes local gain and loss of the stationary 
current density j = as an s-wave ^jJ{r) ‘propagates’ along (the same is 

true for a one-dimensional wave propagating along R [26]). Hereafter 2 ;* means the 
complex-conjugate of z G C. The factorization is also used to transform optical potentials 
with a deep real part Vo into complex potentials with a shallow real part Vo by “removing 
bound states” of the Schrodinger equation associated to U = Vo -|- iVi. The result is a 
class of supersymmetric “phase-equivalent complex potentials” which has been applied 
to study the a + scattering [T8|. The reader can find a thorough discussion of the 
capabilities of the factorization method to intertwine two different complex potentials that 
share the same real spectrum in [T9|. In the same context, the PT-symmetric potential 
V{x) = —{ix)^/2 [50] has been used to generate families of complex potentials which have 
not the property of parity and time reversal invariance but admit a real spectrum ED- 

The organization of the paper is as follows. In Section [2| we use the Ermakov theorem 
to characterize the solutions of the eigenvalue equation Hu = eu, with H and e real, 
in a form that u can be chosen either real or complex-valued. Consistently, these u- 
functions give rise to superpotentials that depend on the parameters of the Ermakov 
equation. In Section [3] we show that new exactly solvable complex potentials with real 
spectra can be generated by using such Ermakov-parameterized superpotentials. It is 
also shown that the left eigenfunctions of the new Hamiltonian H are not equal to its 
right eigenfunctions because this last operator is not self-adjoint. Therefore, although H 
and H are supersymmetric partners, the conventional intertwining relationships are not 
enough to fully describe the properties of the vector spaces involved. It is then necessary to 
extend the operator algebra in order to include the intertwining relationships that connect 
H with H\ the adjoint of H. Using the new (generalized) intertwining relationships we 
show that the eigenvectors of both, H and H\ form a bi-orthogonal system. In Section 0] 
some examples are presented. These include diverse families of new complex periodic VT- 
symmetric potentials as well as new complex-valued Poschl-Teller potentials that have a 
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single bound state fSection 14.ip . We also obtain new complex-valued potentials with the 
set of energies of the harmonic oscillator plus an additional real eigenvalue fSection I4.2p . 
In all cases we use systematically a parametric version of the general solution of the 
Ermakov equation ([T]), the derivation of which can be found in Appendix The paper 
ends with some concluding remarks. 


2 Complex-valued wave functions and corresponding 
superpotentials 

Consider the one-dimensional Hermitian Hamiltonian (in proper units): 

H = -— + Vix), Dom(I/)CM. (3) 

dx^ 

We are interested in solving the eigenvalue equation 

— u''{x) -f- \y{x) — e]u{x) = 0, e e M, (4) 

where f'{x) stands for the derivative of f{x) with respect to x and m is a complex-valued 
function which is not necessarily normalizable in Dom(I/). 


2.1 Complex-valued wave functions 


Let us take 


u{x) = uq exp 


f3{y)dy 


(5) 


with P{x) a comp lex-valued function to be determined and uq an arbitrary complex con¬ 
stant. The introduction of ([5]) in (|1]) leads to the Riccati equation 


- l3'(x) + 0^[x) + {e -V (x)] - 0. 


( 6 ) 


It is well known that the relationship ([5]) is reversed by the logarithmic equation 

^{x) =—^\nu{x), (7) 

so that dZD ‘ linearizes’ the Riccati equation ([HD by transforming it into the Schrodinger 
equation (|1]). As we look for (3 : Dom(I/) —)■ C, it is useful to write (3{x) = (3ji{x) +i(3i{x), 
with Pn and f3i real-valued functions to be determined. Then u is rewritten as m = 
where the amplitude <h(a:) and phase H(a:) are dehned by the real-valued functions 


$( 0 :) = <I)o exp 


f3R{y)dy 


X = 


/^i(y)dy + 


( 8 ) 
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Here $o and Sq are integration constants. A straightforward calculation shows that the 
stationary current density associated with u is given by j = 2^\u\‘^ = —2/3/|mP = vp, 
with V = —2(3 1 the corresponding ‘flux velocity’ and p the probability density [26]. As a 
parameter, v is useful in the study of the resonances belonging to short-range potentials in 
one dimension as well as for distinguishing between bound and scattering states [2^126] . 

The real and imaginary parts of equation ([6D lead to the coupled system 

— -I- -I- e — 1/ = 0, (9) 

+ 2PiPr = 0 . ( 10 ) 

From (1T0|) one immediately gets 

— \\iI3i = 2(3r. ( 11 ) 

This last equation suggest the ansatz 

I3r{x) =--^\na{x), (12) 

ax 

with a{x) a function to be determined. This ansatz provides a dehnite form for the 
solution of (ITTl) in terms of the a-function, 

Piix) = ( 13 ) 

a^{x) 

where A is a constant obtained after integration. To get real-valued functions (3r and /?/, 
it is enough to consider either real or purely imaginary a-functions in Eqs. (IT211 and (1T3|1 . 
besides a real A-constant in Eq. ([T3|). Hereafter we assume that a and A are both real. 

The introduction of (1T2|) and (lT3l) into (|9]) leads to the nonlinear second order differ¬ 
ential equation 

a"{x) = [H(x) - e]a{x) + (14) 

a'^(x) 

which is the Ermakov equation ([1]) for M{x) = V{x) — e and Aq = A^ (for details on 
the possible values of A see the items listed below). The Equation flTT|) was actually dis¬ 
cussed in the context of a nonlinear formulation of time-independent quantum mechanics 
by Reinisch [5T] and its solution was used to construct complex solutions of the time- 
independent Schrodinger equation. From now on the dependence of a on the parameter 
A (equivalently Aq) is implicit. If necessary, we shall make it explicit by writing «(x; A) 
instead of a{x). Using ffT21) and flTdll in (0) yields 

<h(x) = <hoa(a:), S(x) - Sq = -A J a~‘^{y)dy. (15) 


Without loss of generality we take <Foe*“° = mq = 1 to write 


ux{x) 


a{x )exp 


-iX / a {y)dy 


( 16 ) 
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Accordingly, the density p\ = \u\\^ = and the flux velocity v\ = —2/3/ = —2X/a^ asso¬ 
ciated to the comp lex-valued function flT^ are posit ion-dependent. Therefore, equations 
dl]) and ffTOj) can be interpreted as the stationary version of the Madelung’s hydrodynam- 
ical equations [52] , 

Considering the relationship = Aq indicated above, one can write A± = ±\/Ao, so 
that our function (fT6|) is decoupled into the pair of particular integrals of the equation 
(|2]) reported by Ermakov in his theorem j27j : 


u\^{x) = a{x) exp 



=Ff\/Ao / a {y)dy 


(17) 


We have the following classification of the w-functions: 


(I) If Aq = 0 then the pair flTT)) is reduced to the same function ux±=o = and the 
Ermakov equation ffT4)) coincides with the Schrodinger equation (jl]). 

(II) If Ao < 0 then one obtains the purely imaginary numbers A± = v^IAqI, which 
are not included in our approach although they lead to real-valued functions ux^ 
(otherwise, the function f3i in flT^ is not real-valued). 

(Ill) IfAo > 0 then the w-functions are complex-valued. 


On the other hand, we like to stress that (IT^ can be expressed also in the form 


ma = <ho<acos I A / a (//)(i//-I-H q — f$o<asin I A / a {y)dy + 


(18) 


where we have retrieved the arbitrariness of the constants $0 and Hq. A simple inspection 
of fflSj) shows that, up to the factor A in the argument of the sin function, the function 
Im(MA) can be connected with the solution, in amplitude-phase form, of the Schrodinger 
equation discussed in e.g. Iig. In this context, to identify the possible bound state energies 
En of the system, one would impose the Milne quantization condition [121 - H6] : 

N{E) = — / a~‘^{x, E)dx ^ n + I, n = 0,1,..., (19) 

^ J — 00 

to get Im(MA) = 0 as |a:| —?■ -|-oo. However, as our function (lT6l) is complex-valued, the 
latter condition affects also the real part Re(MA), so that 

Ux (—l)"'’''^<ho« as |a:| ^ -fcx). 

Thus, if one would like to describe normalized solutions to the spectral problem of H 
with the M-functions, then the condition flT^ is not enough to cancel them at a: = ±cxo; 
it is necessary also that a{x) —)■ 0 as |a;| —)■ -|-cxo. More precisely, in our approach, the 
cancellation of a at the edges of Dom(l/) is required to get regular w-functions addressed 
to represent physical states. Notice that one can proceed also by cancelling Re(MA) first 
and then determining the conditions to cancel Im(MA). 
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2.2 New real and complex superpotentials 


Assuming that the real function a{x) satisfies ina and A is a real number, the solution 
of the Riccati equation ([6]) is given by 




a'{x) . A 

a{x) ^^a‘^{x)' 


( 20 ) 


In general, this last function is meromorphic in an open set of the complex plane with 
poles defined by the zeros of a in Dom(l/). Note that besides the conventional logarithmic 
term (In a)', our superpotential Px includes an additive nonlinear term regulated by the 
A-parameter that distinguishes the Ermakov equation ([1]) from the linear second order 
differential equation ([2]). Taking into account the items (I)-(III) of the previous section, 
we identify three different classes of superpotentials. 


2.2.1 Conventional superpotentials (A = 0) 

If A = 0 the superpotential (1201) is a real-valued function (3x=q = —(Ina)' producing the 
Darboux transformations already studied in the conventional supersymmetric approaches. 
Thus, if A = 0 the derivation of the transformation functions u in terms of the Ermakov 
equation is unnecessary. 

2.2.2 Real superpotentials (A± = v^|Ao| G C) 

If Ao < 0 then A is purely imaginary and the superpotential (l20l) decouples into the pair 
of real-valued functions 

/3a± =-(Ina)'T (21) 

However, as indicated in the item (II) above, this case is not included in our approach. 

2.2.3 Complex-valued superpotentials (A± = ±\/Ao G M) 

If Ao > 0 then A is a real number and the superpotential (|20|) decouples into the pair of 
complex-valued functions 

/3a± = -(In a)' ± (22) 

Observe that (3x_ = (31^. To avoid confusion with notation used in the previous case, 
hereafter we use /3a± = f3± if A± = ±\/^ G M. Remarkably, the superpotentials (l22i) are 
complex-valued functions even if e and R in ([6]) are real. 

Note that the values of A producing real (complex) superpotentials at the same time 
give rise to real (complex) u-functions. 








3 New complex potentials with real spectrum 


Once we have constructed the complex-valued superpotential that solves the Riccati equa¬ 
tion (|6]) in terms of the Ermakov equation ([1]), we are able to factorize the initial Hamil¬ 
tonian ([3]). Let us consider the hrst order differential operators 

A = --^+l3x{x) and B =-^ + I3x{x), (23) 

ax ax 

with /3 a dehned in fl2(l . As /3 a satishes ([HD, the product of A with B gives 

H = AB + e. (24) 


Hereafter the real parameter A will be different from zero, then the function /3 a will be 
any of the complex-valued functions introduced in Section 12.2.31 The description of the 
transformations associated to the real superpotentials derived in Section 12.2.11 for A = 0 
can be done in terms of the already reported supersymmetric approaches. Then, it is a 
matter of substitution to verify that = B^A^ -|- e = iL, as /3^ is necessarily a solution 
of (|HD- Thus, although A and B are not mutually adjoint they factorize H in the ordering 
established by fIMl) [23] . Now we reverse the order of the factors to get (the dependence 
of i/ on A is implicit): 

H = BA + e = -^ + Vx{x), (25) 

where the new potential Vx{x) is a comp lex-valued function of x that is regulated by the 
A-parameter of the Ermakov equation, 


Vx{x) =V{x)+2/3',ix) 

= H(x)-2-^ln«(x)-*4A^ 
dx^ a-^{x) 


(26) 


Notice that, in contrast with H, the operator H is not self-adjoint since V is complex. 
Indeed = A^B^ + e = —-^ + Vx H. On the other hand, it is straightforward to 
verify that the new potential Vx{x) is associated to the nonlinear (Riccati) equation 

= I3'x{x) + I3l{x) + e. (27) 

The comparison of equations ([HD and (HTh makes clear that, in general, they are not 
invariant under the change (3x —/3a- 

The factorizations (|24D and (|2^ lead automatically to the intertwining relationships 

HB = BH, HA = AH, (28) 


which, as is well known, are useful in constructing the eigenfunctions of H by the action of 
B on the eigenfunctions of H and vice versa, one gets the eigenfunctions of H by acting A 
on the eigenfunctions of H. This point will be dealt with further in the next subsections. 
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3.1 Eigenvalues and eigenfunctions of the new potentials 


li e ^ E and '^e(x) is a normalized solution of the eigenvalue problem HxjjE = E'lpE, 
then 'ipE oc BipE is a solution of HipE = E'lpE that can be normalized for properly chosen 
M-functions (a simple calculation shows that A reverses the action of B). There is only 
an additional eigen-solution '0^ of the spectral problem of H that is not considered in the 
intertwining transformations fl25]) . Moreover, it is associated to the eigenvalue e ^ E. 
Such a function is of the form '0^ = with an arbitrary constant, and satishes the 

condition = 0. It must be clear that the analytical properties of V'e are determined 
by the u-function. For instance, using (I16p and A G M we have 



(29) 


Therefore, the zeros of ux are poles of and vice versa. As A G M, the right hand 
side of (12^ indicates that it is sufficient to analyze the zeros and divergences (if any 
exist) of a in order to get information of the square-integrability of An a-function 
that is free of zeros in Dom(I/) and diverges at the edges of Dom(I/), for example, is 
the ideal candidate to construct 'ijj^ with hnite norm. Considering these remarks, if V'e is 
normalizable in Dom(I/) = Dom(C), then it must be added to the set {iPe}, so that H is 
almost isospectral to H because cr{H) = a{H) U {e}. Otherwise cr{H) = 


3.2 Orthogonality in the new vector spaces 

At this stage we would like to emphasize that although the functions 'i/je can be normalized, 
they are not mutually orthogonal. In this respect, if necessary, one can introduce a 
positive-dehnite scalar product under which the set {iPe} is orthonormal (see for instance 
[53l[5l], ‘atypical models’ in [3], and references quoted therein). Another option arises by 
noticing that the algebraic relationships fl28|) lead to the maps 

B-.n^n and A-.n^H, (30) 

where the Hilbert space Ti = span{'0E} is the state-spaceH associated with El and EL is 
the representation space of H induced by B. As the operator H is not self-adjoint we 
have formally a second pair of intertwining relationships 

HB^ = B^H\ H^A^ = A^H. (31) 

In consequence, there is a pair of additional rules 

A^ and B^ (32) 

^In this section we assume that the spectrum cr{H) of H is purely discrete (j{H) = ad{H), so that the 
inner product {Ee^Ee') = ^ee', with 5ij the Kronecker’s delta, holds for any pair of energy eigenstates 
Ee and ipE'- If scattering states are also present [(Jc{H) E 0); then the Dirac normalization [EetEe') = 
5{E — E') must be considered for E and E' in ac{H). However, some caution is required in order to 
extend the discussion of this section to the states belonging to the continuous spectrum of H [551157] . 
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Here "H is the representation space of induced by Hh Therefore, the action of H, 
and on ‘kef vectors differs from their action on ‘bra’ vectors: 


f B\'iI)e) oc \'iI)e) 

\ A\^e) oc \iIje) 

f {'iIje\B^ oc {^pE\ 

\ {^e\A^ oc (V'eI 

(33) 

1 A^\iIje) oc \iIje) 

J oc (V'eI 

(34) 

1 B^\'tpE) oc \'tpE) 

1 {^^eIB oc {'iI)e\ 

These last right-left action equations, together with the intertwining relationships fl28ll 
and fl3Tll. give rise to two different forms of calculating the matrix elements of the diagonal 
operator H — e. Namely, for E ^ e, the matrix element 

{iI^eHH - e)\'il)E') = 

{E - e){'ijjE\'4^E') 

(36) 

can be expressed as either 



{'i/jeIABI'i/je') oc {i>E\^E') or 

{ijElB^A^li^E') oc {^eIiPe')- 

(36) 


As {'4>e\'4’e') = ^EE', we can use fl35|) and (l36|) to £x the constants of proportionality 
indicated in as follows 



B\i^E) 

^/E^e 


IV^e) — 




(37) 


In this form the vectors |'0 e) and |'0^) are bi-orthonormal 

{'ll>E\^E')fl = ^EE', {^E\'ipE')n = ^EE', {i>E\^E')n = {^E'\'i/Je)h ^ ( 3 ®) 

where the subscripts denote the representation space in which the kets are assumed as 
vectors and the bras as functionals. For instance, the operators BA and A^B"^ are dehned 
to act on "H and "H respectively, then (remember E ^ e): 

{'iIje\BA\^e')u = {^E\A^B^'ijjj^,)^ = - ef\^pE') = {E - e)5EE'- (39) 


The introduction of a bi-orthogonal system like {IV’e), Ii/'e)} dehnes a connection between 
the representation spaces 77 and 77, with the inner product (-I-) in 77 as the mediator. 
Such an structure closes the spectral correlations induced by the algebraic rules (1^ and 
ca. and provides a mechanism to investigate the completeness of the eigenvectors of H 
(see for instance [SH]). Nevertheless, a few precautions are necessary; to name some of the 
related challenges: (i) the physical interpretation of the products fl36|) is not immediate if 
they lead to complex numbers [59]; (ii) the completeness inferred, at hrst sight, from the 
products (l38ll would be not properly justified unless a series of conditions are satisfied [60] . 
see also [61] and (hi) in contrast with the feasibility of hnite-dimensional bi-orthogonal 
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systems [58], the existence and completeness of bi-orthogonal systems in the inhnite- 
dimensional case is not obvious, specially if the continuum spectrum is present [551 - I58] . 
Considering this last point, in particular, the resolution of identity (and thereby the 
quantum averages of observables) could be threatened if vectors l"^) with zero-binorm 
(iplip) = 0 are present [55f[57] . However, such a “self-orthogonality” phenomenon is 
neither accidental nor rare in physics [62], so they must be considered in the approaches 
on the matter. Quite recently, for instance, some of the difficulties involving zero-binorm 
bound states in a continuous spectrum (just as possible divergences in calculating quantum 
averages) have been faced with a proper resolution of identity [55]. In what follows we 
pay attention to a pair of prohles associated with the rules of our approach. 

They correspond to the case when is the complex-conjugate of its dual \^jJE), and 
to the presence of the missing state |'0e) in our model. Major precisions concerning the 
resolution of identity and the calculation of quantum averages of observables will be 
presented elsewhere. 

Using the resolution of identity I = /^ dx\x){x\, the ^-representation of the bi-product 
in Ti is introduced as 




n 


Dom(y) 


/Dom(y) 


1p*E{x)'lpE'{x)dx. 


(40) 


Consistently, {'ipE^'^E')'fi = (V'e'jV'e)^ = ^ee'- In dehnition fl40|) we assume that the 
Darboux transformation (12611 produces potentials V\{x) that are bounded from below, 
just as V{x) is. This can be achieved, for example, if the real part of V\{x) is bounded 
from below and it dominates over the imaginary one such that the ratio \mV\/IieV\ 
remains hnite and sufficiently small at both edges of Dom(U). Complex potentials with 
these properties have been glimpsed as a possible result of the supersymmetric quantum 
design at the time that they were classihed as soft non-Hermitian [56]. As we are going 
to show in the next sections, this is precisely the type of potentials that arises from 
our model because ImlA behaves as a~^ whereas RelA includes an additive term which 
behaves as with a a solution of the Ermakov equation fll4p . Thus, the properly 
chosen a-function allows the calculation of the inner product {gp^,ipE')ii by preserving 
the boundary conditions on Dom(U) C M, as it is indicated in fl4Up . In turn, from (137p 
the ^-representation of \'iIje) and \'4>e) is respectively given by 


'^e{x) = We{x) + P{x)'iI)e{x)] , f^E^x) = We{x) + P*{x)'iI)e{x)] . (41) 

Then, if fjEix) is a real-valued function (as it occurs for the majority of the one-dimensional 
real potentials), we have iP*e{x) = iIje{x). Thus, the wave-function f^Eix) and its dual 
fjEix) are complex-conjugate. In such a case the inner product (jTni) acquires the form 


{'ljj*E,'ljjE') = {'lpE^'^*E')= 'lljE{x)'l/jE'{x)dx= 'ljjE{x)'ljjE,{x)dx = 6EE'- (42) 

JDom{V) JDom{V) 

On the other hand, as indicated in the previous section, the missing state I'l/jf) is in the 
kernel of the intertwining operator A. Therefore, given |'0g) G H with E ^ e, using (IHTD 
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one has 


1 


(43) 




y/E^ 


{i/jE\A\ij,) = 0 . 


Hence, the missing state iV'e) is bi-orthogonal to all the vectors \^jJE^e) ^ Now, let us 
investigate if there is a non-trivial vector \Tp) e H such that it is bi-orthogonal to all the 
vectors ^ From flHTll we get 


(vjIV'e) = {^\B\'iPe) = 0 B^\^p) = 0. (44) 

Therefore = e|^), so that \ip) = is dual to |'0e) and is also bi-orthogonal to all 

the vectors ^ B. Moreover, in contrast with ([39]), which was calculated for E ^ e, 

now we have 

= {'i/j,\A^B^'i/j^)^ = 0. (45) 

This last result certihes that I'lpf,) and are respectively eigenvectors of H and iJl 
with eigenvlue e G M. To investigate the normalization of the missing state notice that 
{x) = 'ipei.x) = c^/ux{x), therefore 


{'ip^ipP) = I 'ipl{x)dx = I ip*J‘{x)dx. (46) 

iDom(y) JDom(V) 

In general, the bi-product (146|) depends on the analytic expression of the a-function that 
dehnes the supersymmetric transformation. In particular, if Dom(I/) C M is symmetric 
with respect to the origin, the a function can be constructed such that fl46|) is a hnite 
real number, as it is shown in the next sections. In this form, the bi-orthonormal system 
{{iV'e), IV’e)}, {iV'e)) |'0 e)}} underlies a viable mathematical structure for our model. 


4 Spectral design of complex-valued potentials 

Let V be an exactly solvable potential such that Dom(I/) C M. That is, V can be 
one-dimensional (Dom(I/) = M) as the free-particle, the modihed Poschl-Teller or the 
harmonic oscillator potentials. This can be also radial (Dom(I/) = M’*') as the Coulomb 
potential, or it can be dehned even in a hnite domain as this occurs for either the inh- 
nite square well, the trigonometric Poschl-Teller or the Scarf potentials. The Darboux 
transformation (1261) and the intertwining relationships (|28|) allow for the construction of 
new exactly-solvable complex potentials V with a spectrum that is either identical to 
or almost the same as the one of the potentials indicated above. As we have seen, the 
conditions for mapping the physical wave functions of V into physical wave functions of 
V depend on the analytical properties of the a-function. Next we show the applicability 
of the method using the free-particle and the harmonic oscillator potentials as point of 
departure. In all cases we use a parametric form of the solution to the Ermakov equation 
that is convenient to simplify calculations, details on the derivation of this expression can 
be found in Appendix 
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4.1 Complex supersymmetric partners of the free particle po¬ 
tential 

For the free particle potential Vfree{x) = 0 one can take 2 ; = and v = as the 

fnndamental integrals of (El) with /c^ = e > 0, and Wronskian W{z,v) = —2ik. Using 
Eqs. flA.4ll and flA. 6 l) of Appendix |A] (with Aq = A^), the parametric family of solutions 
to the corresponding Ermakov equation is given by 


a{x) 


ae-2*^" + ce2‘'=" + 2A/ac + 


4/c^ 


1/2 


(47) 


4.1.1 New complex periodic PT-symmetric potentials 


As fc G M and A G M, we have two immediate (but not exclusive) forms to get a real 
function a. Namely, we can take either a = c = 1/2 or a = —c = In the former case 
one has 

a{x) = [cos{2kx) + 7a], 7a = + (48) 

The superpotential (l20|) and the new potential (l26|) are in this case 


= 

Remark that 7 a > 1 because A 7 ^ 0, then the above expressions are free of singularities. 
The potential Vx{x) has been depicted in Figure [T]for different values of A and k. In turn, 
the square amplitude Pe{x-, A) = iV'eP of the missing state 


ksm{2kx) + iX 
cos( 2 /cx) + 7 a 


Vx{x) = 


4fc^ [1 + 7 a cos(2fcx)] + iAXk sm{2kx) 
{cos{2kx) + 7a)^ 


(49) 


p,(l; A) 


cos( 2 /cx) + 7 a 


(50) 


is a legitimate probability density whenever the integral J^pe{x;X)dx is hnite. Clearly, 
given /c G M, there is no value of A G M such that this last condition is fulhlled. That is, 
using only the physical energies E = k^ > 0 oi the free particle to dehne the n-function 
we obtain a family of complex periodic potentials Va that have no bound states. How¬ 
ever, it is well known that the spectra of periodic potentials are composed of allowed and 
forbidden bands of energy, the former associated with Bloch functions as eigen-solutions 
of the related Schrodinger equation [63]. Then the spectral analysis of the periodic po¬ 
tentials 0491) can be done in terms of the transfer (or monodromy) matrix [63] and will be 
reported elsewhere. On the other hand, the real and imaginary parts of Vx{x) are clearly 
even and odd functions so that this last potential is invariant under the transformations 
i —)■ —i, X — —X. In other words, given A, the complex periodic potential Vx{x) is a VT 
supersymmetric partner of the free particle potential. Compared with other expressions 
of periodic PT-symmetric potentials (not obtained in supersymmetric approaches), our 


14 













result fH9|) could be classified as a generalization of the Scarf I-trigonometric-type poten¬ 
tials reported in e.g. [M]- Interestingly, the periodicity of potentials like the ones defined 
in can be connected with the periodicity of optical crystals that are balanced in their 
gain and loss profiles [55] . 


20 
10 
0 

-10 
-20 
-30 

-10 -5 0 5 10 -10 -5 0 5 10 -10 -5 0 5 10 

(a) k = 0.5, A = 1 (b) A: = 1, A = 1 (c) /c = 1, A = 0.5 

Figure 1: Members of the family of periodic complex potentials V\{x) introduced in (H51) for the indicated 
values of A and k. In all cases the continues (blue) and dashed (red) curves correspond to the real and 
imaginary parts of the potential, respectively. 




The other immediate option a = —c = i/2 produces 

a{x) = [sm{2kx) +, (51) 

with results similar to the ones indicated above. 


4.1.2 New complex Poschl—Teller-like potentials 


To broaden the range of applicability of our approach let us consider unphysical eigenval¬ 
ues {E = k'^ < 0) of the free particle potential. Using k = with k > 0, the periodicity 
of equations flT8|) -fl5T |) is lost. Then, for the results associated with the superpotential 
fH9|) we have the hyperbolic expression^; 


a{x) = [cosh(fi;x) + Ox^^'^ , I3x{x) 


— I sinh(ra) -1- iX 
cosh(fi;x) -b 6x ’ 



(52) 


Note that the a-function is real-valued for A G [—f, f]- The new potential and square 
amplitude of the missing state are given by 


Vx{x) = - 


-b 9x cosh(ft;a;)] -b i2AKsinh(ra) 
(cosh(fi;a:) -b Ox^ 


p,{x] A) = 


cosh(ra) -b Ox 


(53) 


In this case the amplitude is finite along the straight line for all the values of k > 0 
and A G [—I, |]. Indeed, the constant that normalizes such an amplitude is given by 


A 


c, = 


2arctan (1 — 6 *a)] ’ 


with lim Cf = \ l 

A^o V 2 


(54) 


^In this (hyperbolic) case, for simplicity in notation, we omit the tilde over the potential (1531) and the 
wave function of the missing state (1551) . Other functions as a, Px or do not include different notation. 
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Therefore, given A and k, the complex-valued potential Vx{x) dehned in fl53D has the 

2 

single energy E = —^ with eigenfunction 


^pe{x; A) = c^a exp 


|i arctan 


■ K 


(1 — 9x) tanh 



(55) 


The function Vx{x) in (l5^ really represents a A-parameterized family of complex-valued 
potentials that share the same single-point discrete spectrum (t{Hx) = {E = In 

Figure El subhgures (a) and (b), we show two different members of the family Vx{x) that 
share a single discrete energy E = —1/4 (i.e. k = 1). The corresponding wave functions 
'0e(x;A) are included in subhgures (c) and (d). Note the sensitivity of the imaginary 
parts of both the potentials and the wave functions to the variation of A, in all cases their 
absolute value increases with |A|. The real parts, on the other hand, have a less important 
sensitivity to the changes of A. As we are considering A G M, we know that (5- is the 
complex conjugate of /!+ (see Section [2. 2. 31) . The latter property produces V_{x) = V^{x) 
and '0e(a:; A_) = 'il}*{x] X+), so that the imaginary parts of V_{x) and '0e(a:;A_) are the 
negative version of the ones depicted in Figure El Indeed, the former correspond to the 
mirror image of the latter and vice versa. 



Figure 2: (a-b) Members of the family of complex-valued potentials Vx(x) defined in (1551) for k = 1, 
A = 0.2 (black, dashed) and A = 0.45 (red, continuous). The dotted horizontal line, depicted in the 
real part of Va, refers to the single point discrete spectrum {E = —1/4} that is shared by all the family 
members, (c-d) The corresponding wave functions il)e{x\X). 


• "PT-symmetry 

The symmetry indicated above is clear for since the argument of the exponential in 
(155|) is pure imaginary and the function arctan(a:) is odd. Therefore, the change A —A 
is equivalent to changing i —i. Concerning the potential Vx{x) introduced in (15^ . 
let us emphasize that the change x —x reverts the i-transformation because a and 
cosh(a;) are even functions while sinh(a;) is odd. That is, the family of potentials reported 
in fl5^ is PT-symmetric as the transformation x —x, i —)■ —i, leaves it invariant. 
The same holds for Note however that the superpotential fix becomes —/dy after such 
a transformation. We want to stress that our potentials (1551) could be classihed as a 
generalization of the Scarf I-hyperbolic-type potentials reported in e.g. [51]. A similar 
classihcation holds for the potentials that are derived from the hyperbolic version of the 
superpotential (l5TD . 

• Supersymmetric Poschl—Teller profile 

The explicit analytical form of the new potentials fl53|) reveals a prohle that is well 
known in the supersymmetric approaches associated to the free particle potential miEsi 
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EZ]. As a matter of fact, the modified Poschl-Teller potential Vpt{x) = cosh ‘^(nx) is 
one of the regular supersymmetric partners of Vfree{x) = 0 [ 66 ]. This last result gives rise 
to the identihcation of the delta well potential Vs{x) = —jo6{x), with 70 a properly chosen 
constant, as another supersymmetric partner of Vfree [SZ|- Here, the family of complex 
potentials is also of the Poschl-Teller form. To conhrm such statement remember 
that A = 0 must reduce our results to the conventional ones fSection l 2 . 2 .ip . This last 
is clear since a{x; \ = 0 ) = [cosh(fi;a:) + 1 ]^/^ = \/ 2 cosh(^), which leads to the regular 
superpotential (3x=o{x) = —|tanh(^) reported in [ 66 ]. Hence, the entire family Vx{x) is 
in the class of regular supersymmetric partners of the free particle potential and has a 
Poschl-Teller prohle (compare with e.g. [T^EU] ). 

To give a particular example let us cancel the number 6x dehned in ([52]). Then, using 
A = I in fl53|) . we arrive at the function 


Vk{x) 


.K^sinh(fi:a:) 
cosh^(«:a:) cosh^(ra) 


[1 + isinh(ra)] Vpt{x). 


(56) 


This last expression includes the modihed Poschl-Teller potential Vpt{x) as its real part 

2 

and has a single bound energy e = — ^, with normalized wave function 




K 


7rcosh(ra) 


1/2 


exp 


iarctan 


tanh 



(57) 


The potential Vk{x) and the wave function 'ipeix) have been depicted for different values 
of K in Figure [3l 


12 

0 


-25 

(a) Fk, k = 5 ("pj ^ _ 5 (c) Vk, K = 2 (d) V'ej ^ 5 






Figure 3: Real (blue, continuous) and imaginary (red, dashed) parts of the complex-valued Poschl- 

Teller potential Vk{x) and the wave function for (a-b) k = 5 and (c-d) k = 1/2. The single discrete 
^ 2 

energy e = (dotted, black) of Vk(x) is less negative than the single energy Eq 1.5e (black) of the 

2 

modified Poschl-Teller potential Vpt(x). The imaginary part of the potential is shadowed to emphasize 
its gain and loss profiles. 


• Optical model 

Using the optical model discussed in the introduction one can associate the imaginary 

part VAx) = sinh(fi;a:)l/pT(a^) of the potential Vk(x) with the absorption and emission of 

2 

probability waves in (0,+cxo) and (—cxo,0) respectively (see Figure [3|). Notice that the 
rates of these processes coincide at x and —x respectively since the function Vi is odd, so 
that the total probability is conserved. In Figure [3] we have shadowed (in light orange) 
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the area from Vj to the axis to emphasize the zones of Dom(VK) where this function 

2 

works as a sink {Vi < 0) or as a source (V/ > 0) of probability waves. In this context, 

the complex potential Vk corresponds to a modihed Poschl-Teller system VpT embedded 

2 

in an environment Vj which preserves the number of particles because the source and 
sink operate at the same rate in (—cxo, 0) and (0, oo) respectively. For scattering energies 
the function Vpt{x) represents the average potential to which an incoming particle is 
subjected and Vi{x) regulates the ‘capturing’ and ‘releasing’ of such a projectile by the 
scatterer in a vicinity of the maximum interaction. As Vk{x) is "PT-invariant, the analysis 

of the scattering properties of this potential could be completed by following e.g. Ref. [68]. 

On the other hand, it is well known that Vpt{x) has a single bound state of energy 
Eq = — 1)^ ~ 1.5e [67|. Then, as the single energy of Vk{x) is at e = —the 

effect of the environment Vj on the Poschl-Teller system Vpp is the shifting of the single 

2 

bound energy from Eq 1.5e to the less negative value e = — x- 


4.2 Complex potentials with energies of the harmonic oscillator 

As the fundamental solutions of the harmonic oscillator potential Vosc{.x) = we use 


z{x] e) = lEi 


1-e 1 


4 ’ 2 


,x 


e v(x; e) = XiEi 


3-e 3 


4 ’2 




,-xV2 


(58) 


with iFi[a, c; z\ standing for the confluent hypergeometric function [69117T] and Wronskian 
W{z,v) = —The introduction of these functions in the Equation (IA.4|) of Appendix lAl 
gives rise to an a-function that is parameterized by four number^: a, b, c and e. This 
means we have at our disposal a vast source from which diverse families of complex 
potentials can be constructed with the energy spectrum of the harmonic oscillator. Indeed, 
as the /^A-function fl20|) is always complex and its real and imaginary parts include the 
function a in the denominator, we have to choose the parameters a, b, and c such that 
a is free of zeros in Dom(VA) = for a given e. This last prevents f3x, V\ and V’e of 
having singularities and opens the possibility of getting a well dehned transformation of 
the initial wave functions into the new ones ipp ^ As indicated in Section IXTl we will 

hnd that V\ is strictly isospectral to Vosc whenever is not normalizable and (3x is free 
of singularities. However, our interest here is to manipulate the spectrum of the initial 
potential in order to design the spectrum of the new one. Henceforth, we look for a set 
of parameters {a, b, c, e} leading to regular complex superpotentials (3x and normalizable 
missing states ipe- As an example let us consider the case e = — 1, we have 


z{x) = z{x^ —1) = f (x, —1) = Erf(x) (59) 

"^Note that according to Eqs. (IA.5I) - (IA.7I) of Appendix the parameters a, b, and c are respectively 
equivalent to the constants cq and ci that define the general solution of the Schrodinger equation (^1). 
and to the parameter Aq that characterizes the Ermakov equation O- These last parameters are used in 
Eqs. (lA.ll) and (IA.2I) to define the a-function in Appendix 
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with Erf(x) the error function [^. Therefore, the a-function flA.4D can now be read as 
follows 

a{x) = \a Erf^(x) + b Erf(x) + c] , (60) 

where unnecessary factors have been absorbed into a and b. 


4.2.1 Unbroken snpersymmetry 


Without loss of generality, let us assume that the parameters a, b and c in (160]) are all 
non-negative. As Erf(a:) G [—1,1] we realize that the a-function is real and free of zeros 
if c > 1^. Moreover, a as |a:| —)■ -|-cxo, so that ^ vanishes at the edges of 

Dom(VA) = R- Thus, for e = —1 we hnd that any set of non-negative parameters 
{a,b,c > 1^} produces regular superpotentials f3x and normalizable functions 'ipt- 
other words, the spectrum of the family of potentials 


Vosc{x; A) = Vosc{x) + 2j3'y^{x) 


= X 


2-2 


A 

dx 


b -|- 2aErf(a:) — zyOrA 


(61) 


is the set of energies aiVosc) = {e = U cr(V),sc) = {“1; 1)3,...}. In the usual nomen¬ 
clature one says that each member of the family Vosc{x', A) is a supersymmetric partner of 
the harmonic oscillator potential Vosc{x)- Moreover, as the potentials belonging to such 
a family are almost isospectral to the harmonic oscillator, one also says that every pair 
(Vosc, lose) corresponds to a system with unbroken supersymmetry [Tl fTOlITT] . 





(a) Vosc{x; A) (b) Vosc{x) (c) A) 

Figure 4: (a-b) The supersymmetric pair of oscillator potentials {Vosc, Use) with unbroken supersym¬ 
metry generated with the parameters a = b = and c = 1 (equivalently, A = The spectrum of 

the complex-valued potential Use is real and includes the energy e = — 1 that is missing in the spectrum 
of the oscillator potential Use- The imaginary part of Use is in (dashed) red and shadowed to mark the 
zones in which it works as either a source or a sink. The TT-transformation of Use produces the figure 
(c), so that this potential is not TT-invariant. The potential U*se(~^!A) is, however, a member of a 
family of potentials for which b —>■ —b, so that it is also a supersymmetric partner of Use- 

Ill Figure S] (a-b) we show the behaviour of the supersymmetric pair (Uso Use) for 
specihe values of the parameters a, b and c. These potentials share the same energy-levels 
E = 1,3,5,..., but the lowest one e = —1 is included only in the spectrum of Use- 
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Following Section |3l one can depart from the spectrum of Vosc to obtain the one of Vosc 
by the action of the operator B (just as we have proceeded in the above section) and 
vice versa, the action of the operator A on the wave functions of Vosc reverts the previous 
operation. The first four wave functions of the new potential Vosc are depicted in Figure [5] 
together with the oscillator ones. The missing state iJeix), illustrated in Figure E] (a), has 
no supersymmetric partner in the set of eigenfunctions of Vosc{x) because the spectrum 
of this last potential does not include the eigenvalue e = — 1. Notice that in all the cases, 
the sign of the imaginary part of the potential and the wave functions changes around 
the minimum of the real part of the potential. 

We want to remark that the ‘complex oscillator’ that is usually reported in the lit¬ 
erature is a variation of the conventional oscillator in which the frequency is allowed 
to be a complex number (see e.g. [MZH and references quoted therein). As a conse¬ 
quence, the eigenvalues of the related non-Hermitian Hamiltonian are complex numbers, 
so that the spectrum is highly unstable under small perturbation^ [73]. In contrast with 
the complex-frequency oscillator indicated above, our complex oscillators fl6T]) are free of 
spectral difficulties because their frequencies and eigen-energies are purely real. 



Figure 5: The first four wave functions of the complex-valued potential Vosdx] A) shown in Figure|3](a). 
In all cases the real part is in blue and the imaginary part in (dashed) red. Figures (b-d) include (in dotted, 
black) the corresponding oscillator wave functions. The wave function in (a) has no supersymmetric 
partner in the set of eigenfunctions of Vosdx) because the energy e = — 1 is not included in the oscillator 
spectrum. The wave functions associated to the FT-transformed potential of Figure|4](c) are immediately 
obtained from the ones shown here by using the same transformation. 


4.2.2 Broken PT-symmetry 

For simplicity, in the previous section we have used a set of non-negative parameters 
{a, 6, c > 1^} to get regular Darboux transformations of the oscillator spectral problem. 
Of course, this last set is not unique but serves to facilitate the identification of the new 
families of complex potentials with an exact solution to their spectral problem. Remark¬ 
ably, such a set of parameters can be extended by the symmetries involved. For instance, 
the change 6 —)■ —6 preserves the analytical properties of the a-function and gives rise to 

^Despite this difficulty, a Lie-admissible theory has been applied in the study of the density matrices 
that represent the quantum states of such a complex oscillator m and numerical approximations have 
been developed to investigate the related resonances m- Moreover, the corresponding supersymmetric 
partners have been used to obtain solutions to the Painleve IV equation m- 
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the potential depicted in Figure |4] (c). Thus, the change b ^ —b corresponds to the VT- 
transformation x —)■ —x, i —)■ —i, so that it produces Vosc{x; X) —)■ V*^^{—x;X). This last 
remark means that the new set of parameters {a, —b, c > |^}, with a and b non-negative 
numbers, represents another family of complex-valued potentials with a spectrum that 
is almost the same as the one of the harmonic oscillator. A member of this last family 
is precisely the potential in Figure 0] (c). Considering the results of the Poschl-Teller 
families discussed in Section 14.11 we realize that the hnding of symmetries in the set of 
parameters obtained from the Ermakov equation to dehne the Darboux transformation is 
quite natural. However, in contrast with the periodic and Poschl-Teller cases discussed 
in the previous sections, here we have found a mechanism to produce complex potentials 
with real spectrum that are not invariant under PT-transformations. 


5 Concluding remarks 

We have presented a new way of constructing complex superpotentials in supersymmet¬ 
ric quantum mechanics. The model is based on the relationship between the nonlinear 
equation ([1]) introduced by Ermakov and the linear second order differential equation (|2]). 
The new aspect here is mainly represented by the transformation function ux{x), used to 
generate the superpotential = —(Inu)', because this is parameterized by the solutions 
of the Ermakov equation. In this form, the superpotential I3x{x) is for sure a complex¬ 
valued function, no matter if the initial potential V (x) and the factorization constant e 
are real in the Riccati equation (|6]) from which the supersymmetric transformation orig¬ 
inates. As usual in the supersymmetric approaches, the method presented here includes 
both almost- and strictly-isospectral partner potentials. In all cases our model ensures 
that the new complex potentials Vx{x) are bounded from below in the context of the soft 
non-Hermiticity studied in |nS] , and allows the construction of a bi-orthonormal system 
formed of the eigenvectors of Vx{x) and the ones of its Hermitian conjugate. Moreover, 
in the almost-isospectral cases, the appropriate solution of the Ermakov equation gives 
rise to a vector state '0e(x) that is bi-orthogonal to the other eigenvectors of the 

complex potential Vx{x). The function '^^{x) represents the ground state of Vx{x) and 
belongs to the real energy e, which is missing in the spectrum of the initial potential 
V{x). Further properties of this pretty special ‘missing state’ and the completeness of the 
bi-orthonormal system associated with Vx{x) are to be investigated. 

We have illustrated the method by constructing new supersymmetric partners of the 
free particle potential that are either complex PT-symmetric periodic potentials or com¬ 
plex PT-symmetric regular potentials with a single discrete energy. The latter are of 
Poschl-Teller form, as expected for regular supersymmetric partners of Vfree{x) = 0. We 
have constructed also a family of complex-valued potentials Vosc{x) with the entire spec¬ 
trum of the harmonic oscillator potential {E = 1,3,5,...), plus an additional energy at 
e = — 1. All the members of this family are supersymmetric partners of the oscillator 
potential in such a way that each pair (Wsc, Wsc) forms a supersymmetric system with 
unbroken supersymmetry. Remarkably, the new oscillator complex potentials Wsc have 
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real frequencies and are not invariant under PT-transformations. 

The procedure indicated above can be repeated at will, in order to produce new families 
of potentials. Moreover, it can be applied in the study of other potentials with domain 
of dehnition Dom(l/) C M, including the ones which have been already analyzed with the 
usual supersymmetric approaches. 


A Parametric solutions of the Ermakov equation 


According to Ermakov [27], the equation ([T]) has the solution 


a{x) = z{x) 


Ao 

- \- Cl 

Cl 


z {y) dy + 



(A.l) 


where Ci and Cq are integration constants while z(x) is a solution of the homogeneous 
equation ([2|). The structure of (lA.lD allows to identify the function 

nx 


v{x) = z{x)q{x), with q{x) := Wq 


z ‘^(y)dy 


(A.2) 


as the conventional expression for a second solution of (Ej) providing z{x) as given 
It is assumed that the Wronskian W{z,v) = zv' — z'v is a constant, W{z,v) = tco, so 
that z{x) and v{x) are linearly independent. In consequence, a in flA.l|) is reduced to the 
general solution of ([2|) if Aq = 0, as Eqs. ([T]) and ([2]) coincide for this value of Aq. In terms 
of the fundamental pair 2 ; and n, the solution flA.lIl reads as follows 


a{x) = z{x) 


Cl^2t^\ , O'"0 I -^0 + Cq 


L'«^0 


q\x) + 2^q{x) + 


Wo 


Cl 


1/2 


*^1 1 I o ^0 ^ 1 ^ 1 I / -^0 + Cq 


L '«^0 


v^{x) + 2— v{x)z{x) + 


Wo 


Cl 


1 1/2 


z‘^{x) 


(A.3) 


The discriminant of the square root in the hrst line of this last equation is equal to —4^, 
therefore the solutions of ([ID are elements of the parameterized family 


a{x) = [av'^{x) + hv{x)z{x) + cz^{x)\ 


1/2 


(A.4) 


where only two of the three parameters a, b and c are independent of each other. For 


Cl 

a = c = 

wi 


-^0 + Cq 
Cl 


(A.5) 


it follows, with 


W{z,v) = wo, 6^ — 4ac = —4— 

wi 


(A.6) 
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that 



(A.7) 


where Aq is still involved as a third parameter. The above expressions are consistent with 
the ones reported in e.g. [221130] (see also |12H10])- Remark that any selection of the 
parameters a, b and c prodncing the linear combination a = z + rjv, with 77 an arbitrary 
nnmber, leads to Aq = 0. According to onr discussion in the previous sections we have the 
constraint Aq > 0, so that a will be not the general solution of the Schrodinger equation 
(l 2 |) for the parameters a, b and c that we are considering in this work. 

Another point of interest is to note that the expression 



(A. 8 ) 


with the constant jo introduced for the sake of units, is ^-independent. That is, J' = 0. 
In this context, J resembles the Ermakov dynamical invariant I. While / is a constant 
of motion in the dynamical time-dependent case, the expression J is here just a constant 
with respect to x. If u is of the form ffT 6 |) . then such a constant is equal to zero. Therefore 



with A± = ±\/Ao, see Section 1221 This last result makes clear that u and a are equivalent 


if Ao = 0, as expected. 
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